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It has been recently realized that in many random 
physical processes the conceptions of Gaussian distri- 
bution and Fickian diffusion are invalid. Many such 
processes can be described by models of sub- or super- 
diffusion, where the displacement moments of the corre- 
sponding random walk grow slower or faster than for nor- 
mal diffusion, respectively. A typical example of super- 
diffusion is the enhanced transport in fluids, predicted 
for flows with velocity correlation functions slowly de- 
caying in space or time [1]. A specific type of super- 
diffusion, the Levy flight, has been reported in observa- 
tions of transport in two-dimensional rotating flows and 
in a freely decaying two-dimensional turbulent flow [5]. 
Other examples of super-diffusive transport include wave 
turbulence, non-local transport in plasma, transport in 
porous media, surfactant diffusion along polymer chains, 
cosmic rays propagation and motion of animals [3]. A 
widely used description of super-diffusive transport relies 
on the continuous time random walk model with a power 
law asymptotics of the particle jump length distribution, 
leading in the macroscopic limit to a diffusion equation 
with the Laplacian replaced by its fractional power [4^ . 

An important problem is the influence of super- 
diffusion on processes with chemical reactions [5j [B] . Nor- 
mal reaction- diffusion systems exhibit different types of 
instability [7]. Profound understanding of pattern for- 
mation and spatio-temporal chaos in these systems was 
achieved through generic equations valid near the insta- 
bility threshold, such as complex Ginzburg-Landau ^ 
and Kuramoto-Sivashinsky equations [T]. The evolution 
of instabilities in reaction-diffusion systems can be ac- 
companied by advection of components. For instance 
stirring, which changes the effective diffusion properties 
of species, is one of the means to control dynamical 
regimes generated by instabilities in reaction- diffusion 
systems [S] • Thus one can expect that in some cases flows 
can give rise to an enhanced diffusion of reagents. While 
studies of instabilities in systems with sub-diffusion have 
started (see [10] and references therein), super-diffusive 
reaction-diffusion systems are still unexplored, with the 



exception of the front propagation phenomenon, which 
is strongly influenced by fluctuations [B]. In this letter, 
weakly non-linear dynamics of a reaction-diffusion sys- 
tem characterized by Levy flights near a long wave bifur- 
cation point is investigated. 

Consider a two-component reaction-diffusion system 
in the general case of distinct anomaly exponents: 
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where nj,dj and fj are the species concentrations, dif- 
fusion coefficients and general kinetic functions, corre- 
spondingly. The fractional operator of order 1 < 7 < 2 
is defined as [11] 
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The equivalent definition in Fourier space allows for a 
simple generalization of the operator to higher spatial 
dimensions: 2)j^|e"*'^ = — |q|'''e''i '' . Suppose that there 
exists a homogeneous steady state no satisfying f(no) = 
0. A vanishing trace of the sensitivity matrix, (Vf )jfe = 
dfj/drik, j,k e {1,2}, leads to Hopf bifurcation at the 
long wave limit q = 0. Take e <C 1 and < ~ 0(1) so 
that tr Vf|no = e^ii and invoke a multiple scales analysis 
with n{x,t) = N{^,to,t2,...;e), f = S{e)x, tj eH, j = 
0,2,... and 
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For normal diffusion ("f — 2) S = e, Sj = e-' and 
a sequence of problems at successive orders dj is ob- 
tained. The solution at order 61 is of the form Ni — 
A{£^, t2, ■ ■ .)e*'^°*°vi -I- c.c, where vi is an eigenvector of 
the linearized problem and ujq is the Hopf bifurcation 
frequency. Neglecting the phenomena evolving on time 
scales longer than t — t2, the equation for the ampli- 
tude A ensues as a solvability condition at order ^3. For 
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an anomalous system the scaling property of the frac- 
tional operator, 2)j^|y(a;) = S^D'^^^y{^/6), determines the 
scale of the slow spatial variable, S. Namely, in a more 
common case with 71 = 72 = 7, the scale is (5 = e'^/^ 
and Sj = . The amplitude equation has the form of a 
fractional complex Ginzhurg- Landau (FCGL) equation: 
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(in rcscaled form). This equation was formerly derived 
in |12j in the problem of nonlinear oscillators' dynamics 
with long-range interactions. The parameters a and /3 co- 
incide with those of a normal reaction - diffusion system 
but the Laplacian is replaced by the fractional operator. 
If 71 ^ 72 , the super-diffusion term with the larger index 
is negligible in the long-wave region and 5j = for j <2> 
only. Higher-order powers are fractional and depend on 
the ratio of the anomalous exponents. Then the appro- 
priate scaling is 5 = e^/^ with 7 = min{7i,72}, and the 
expressions for a and (3 are obtained by taking (i2 = if 
71 < 72 and di = if 71 > 72. 

The integro-differential equation Q retains the basic 
symmetries of a normal complex Ginzburg-Landau equa- 
tion (with respect to time and space translations and 
the phase change A i— s- Aexp(iz9)). It is interesting that 
its solutions in the form A(^, r) = B(C) e*(«^-"^), with 
0', w G M, have a symmetry similar to that found by Ha- 
gan [13j . If a solution of this type is known for a pair 
(q!,/3), the solution for a new pair (a',/3') located on one 
of the curves {a — (3)/{l + afi) = const can be found by 
the transformation B = aB' , ^ = b£^', where 



The constant c is chosen so that T converges. Then 
dT/dr = -2(1 + P^)-^JZ^ \dA/dtfd^ < 0, and the 
system relaxes to a certain stable "stationary" solutions 
(the original variable A oscillates with the frequency /3). 
Now consider the traveling wave solutions of Q, 

A, = v/r^M^e^^^«-"^\ - /? - (/3 - a)\q\\ (8) 

A small perturbation a(^, r) about Aq comprises longi- 
tudinal and transverse waves of the form 

a = Ag+fc(r)e^(«+'^«)«+^'=''" + A,_fe(T)e'(«-'=«)«-*'='''', (9) 

with k^,kjj being the respective wave numbers. The so- 
lution (|8| is neutrally stable with respect to disturbances 
— = 0. Further insight into long perturbations re- 
veals that for 0{k^/q) ~ 0{kr,/q) ^ o(l) to leading order 
the growth rate of Aq±k ~ exp(Ar) satisfies 
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Therefore all solutions ([S]) are unstable if 1 -f a(3 < 0, 
i.e. the Benjamin-Feir criterion for a normal CGLE is 
recovered. However, ifl + a/3 > 0, a 7-dependent set 
of unstable wave vectors exists, generalizing the Eckhaus 
instability criterion; 



km I < kl < 1, 
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and the new wavenumber and frequency are q' — bq, oj' = 
a'-6^(l + a'')(a-w)/(l + a2). 

In the special case a = /3 eq. Q, after the phase 
shift A ^ v4exp(— j/Jr), eq.Q, like a normal complex 
Ginzburg-Landau equation [14j , can be written in a vari- 
ational form, 
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where T — U{^,T)d^, and 
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No new instability criteria emerge in the opposite limit 
q ^ fc^, fc^ <C 1. In particular, the spatially-homogeneous 
oscillation Aq — exp(— i/3t) is unstable in the same region 
1 + aP < with respect to disturbances whose wave 
numbers k satisfy 



< |A:r < 



2(1 + a^) 
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The evolution of perturbations near the domain bound- 
ary is expected to be described by an analogue of the 
Kuramoto-Sivashinsky equation [T]. Define 1 -I- a(3 — 
-e, < e < 1, rewrite with x = e^^^S, and T2 = e^r, 
take A — exp (— i/3r2/e^r(x, T2) exp [«(p(x,T2)], where 
r = l + Eilie-'''j(x>'r2), ^ = Ejli '/'j(x,T2), and ex- 
pand exp (±i(p) to obtain the phase diffusion equation at 
order 0(e"^) that, after rescaling, has the following form 
(notations for the rescaled space and time variables are 
the same): 
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The operator (D? i)^ is defined in Fourier space by 



;\x\' 

g|27g*<2x and cannot be simply related to 



3 











1000 




} 






1000 










J 1,02 
1 1.01 



lose, 



Figure 1: Spatio-temporal dynamics of solutions of eq.( |13[ l for 
7 = 2.0 (upper), 7 = 1.7 (middle), and 7 = 1.6 (lower). 



Figure 3: Spatio-temporal dynamics of solutions of eq.Q for 
(a) a = -1, /3 = 1.33, 7 = 2.0; (b) a = -1, /3 = 1.2, 7 = 1.6; 
(c) contiimation of (b) for the same parameter values. 
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Figure 2: (a) Spatio-temporal dynamics of solutions of eq.( 13 1 
for 7 = 1.5. (b) Solutions of (131 at successive moments of 
time. 



the operator 



as the order 27 exceeds the defini- 



tion range in (pk. Eq.(13l is the fractional Kuramoto- 
Sivashinsky equation. 

Fig|l]shows spatio-temporal dynamics of the numerical 
solutions of eq.(13) obtained by means of a pseudospec- 



tral code, using periodic boundary conditions and start- 
ing from small-amplitude random data. The upper figure 
shows the dynamics for 7 = 2 corresponding to a normal 
KS equation: this is a well-known spatio-temporal chaos 
exhibiting merging and splitting of " cellular" structures 
Pj. The middle figure corresponds to 7 = 1.7. One can 
see that along with the chaotic dynamics of " cells" large- 
amplitude traveling "shocks" develop that emit cells still 
displaying chaotic dynamics. With further decrease of 7, 
the shocks appear more frequently, propagate faster and 
their amplitude grows, see the lower figure corresponding 
to 7 = 1.6. When 7 decreases below a certain threshold 
that depends on the domain length, a single traveling 
shock is formed in the whole domain. An example of 
such a shock is shown in Figj2] Here the shock is trav- 



eling with a constant speed (Fig|2^) while its "wings" 
exhibit spatio-temporally chaotic modulations, (Fig|2]3). 
Decreasing 7 results in the increase of the shock ampli- 
tude and after certain critical 7 the shock starts accel- 
erating with its amplitude growing exponentially. The 
shock amplitude grows with the size of the computa- 
tional domain. An asymptotic analysis carried out for 
large-amplitude solutions of ([l3]) shows that the solution 
is of the form (/) — a(r)/(^ — C('''))i where / is an odd pe- 
riodic function and a(r) grows exponentially (despite the 
problem non-linearity) with a certain dependence on the 
domain size and the velocity dS^{T)/dT proportional to 
a(r). The numerical simulations confirm the asymptotic 
analysis. 

Next, numerical simulations of the FCGL equation 
([4| in ID have been performed for the phase turbulence 
regime. Fig|3^ shows spatio-temporal dynamics typical 
of the normal CGL equation, starting from the Benjamin- 
Feir-unstable, spatially-homogeneous oscillations: it is 
well described by the normal KS equation (see FigjlJ^). 
Figj3|D shows the similar dynamics of eq.(|4]) for 7 = 1.6. 
One can see that, after some period of phase turbulence, 
accelerating shocks form that trigger the transition to 
defect turbulence shown in FigjSfc. The formation of the 
accelerating shocks seen in Figs|3]3,c is consistent with 
the formation of shocks in the FKS equation discussed 
above. 

Fig|4]shows the spatio-temporal dynamics of numerical 
solutions of eq.(|4]) corresponding to the defect turbulence 
regime emerging from the Benjamin-Feir- unstable wave 
^ with q = 0.5 for 7 = 2.0 and 7 = 1.1 (Figsjl] (a) 
and (b), respectively). One can see that in the anoma- 
lous case the defect turbulence has a stronger phase- 
turbulence component and does not consist of propagat- 
ing holes. 

Finally, numerical simulations of FCGL eq.([4]) in 2D 
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Figure 4: Spatio-temporal dynamics of \A\ - solution of eq.Q 
for a = 1, (3^ -1.3 and (a) 7 = 2.0; (b) 7 = 1.1. 

have been performed for the parameter values corre- 
sponding to the formation of spiral waves in the normal 
CGL equation. Periodic boundary conditions and small- 
amplitude random initial data were used. The results 
are shown in Fig|5] One can see that for 7 close to 2 (see 
FigsjS^jd) the formation of a spiral wave is still observed. 
With the decrease of 7 the spiral-wave regime is replaced 
by a defect chaos, however, remnants of the spiral waves 
still can be seen (Figsj5}D,e), with each "spiral" occupying 
a small domain with the domain walls partially melted. 
Further decrease of 7 results in the decrease of the num- 
ber of defects, the domain walls are almost completely 
melted (Figsj5|;,f), and the local wavenumber created by 
each defect decreases. A "phase diagram" of the new 
dynamical states described above in the parameter space 
will be presented elsewhere. 

In conclusion, we have derived fractional Ginzburg- 
Landau and Kuramoto-Sivashinsky equations that de- 
scribe weakly non-linear dynamics of a super-diffusive 
reaction-diffusion system, characterized by Levy flights. 
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Figure 5: Snapshots of solutions of eq.Q for a = 1.5, (3 = 
-0.6 and 7 = 1.9 (a),(d); 7 = 1.8 (b),(e); 7 = 1.05 (c),(f); 
upper figures - Re (A), lower ones - |^|. 



and studied some of their solutions analytically and nu- 
merically. We have shown that super-diffusion can lead 
to a transition from phase- to defect turbulence and to 
destruction of spiral waves. We note that investigating 
the effects of fluctuations on non-linear dynamics of in- 
stabilities in a super-diffusive reaction-diffusion system 
would be of interest since, as shown in [B], the fluctu- 
ations can have a profound influence on the non-linear 
behavior in such systems. However, this topic is beyond 
the scope of the present paper. 
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